Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 



(MN MfeK style file vl.4) 



The properties of field elliptical galaxies at intermediate 
redshift. II: photometry and spectroscopy of an HST 
selected sample 

T. Treu, 1 ' 2 ' 4 M. Stiavelli, 2 P. M0ller, 3 S. Casertano 2 and G. Bertin 1 

1 Scuola Normale Superiore, P.za dei Cavalieri 7, 1-56126, Pisa, Italy 

2 Space Telescope Science Institute, 3700 San Martin Dr., Baltimore, MD 21218, USA 

3 ESO, Karl-Schwarzschild Str. 2, D85748, Garching bei Miinchen, Germany 

4 present address: California Institute of Technology, Astronomy 105-24, Pasadena, CA 91125, USA; e-mail tt@astro.caltech.edu 



1 February 2008 



ABSTRACT 

A sample of field early-type galaxies (E/S0) at intermediate redshift (z ~ 0.1 — 0.6) is 
selected, based on morphology and colours from HST-WFPC2 parallel images. Pho- 
tometric structural parameters (effective radius Re and effective surface brightness 
SBc) are derived through the F606W and F814W filters, using luminosity profile fit- 
ting and two-dimensional fitting techniques. The combined parameter that enters the 
Fundamental Plane (logi?e ~ PSBe, with (3 « 0.32) is shown to suffer from signif- 
icantly smaller uncertainties (r.m.s. 0.03) than the individual structural parameters 
(e.g. ~ 15 per cent r.m.s. on the effective radius). 

High signal-to-noise intermediate resolution spectra, taken at the ESO-3.6m, yield 
redshifts for 35 galaxies and central velocity dispersions for 22 galaxies. Central veloc- 
ity dispersions are derived using a library of stellar templates covering a wide range of 
spectral types, in order to study the effects of templates mismatches. The average ran- 
dom error on the central velocity dispersion is found to be 8 per cent and the average 
systematic error due to template mismatch is found to be 5 per cent. The errors on 
velocity dispersion measurement and the effects of template mismatches are studied 
by means of extensive Montecarlo simulations. In addition, we investigate whether 
the determination of the velocity dispersion is sensitive to the spectral range used, 
finding that the value of velocity dispersion is unchanged when the spectral regions 
that include the absorption features Ca HK and NaD are masked out during the fit. 

Key words: 

galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: photometry — 
galaxies: kinematics and dynamics — galaxies: fundamental parameters — galaxies: for- 
mation 



1 INTRODUCTION 

This paper is the second of a series devoted to the investiga- 
tion of the internal and structural properties of field early- 
type galaxies (E/S0) at intermediate redshift (0.1 < z < 1), 
and of the properties of their stellar populations. In partic- 
ular, the aim of this series of papers is to investigate how 
empirical scaling laws such as the Fundamental Plane (FP; 
Djorgovski & Davis 1987; Dressier et al. 1987) and the Ko- 
rmendy (1977) relation evolve with redshift. These studies 
are complemented with the measurement of metallicity and 
age of stellar populations by means of the Lick/IDS absorp- 
tion lines indeces (Trager et al. 1998). 

For this project we have collected high signal-to-noise 



(S/N) spectra at intermediate resolution (resolving power 
A/AA ~ 1000, where A A is the full width half maximum res- 
olution) for a well-defined sample of field early-type galaxies 
in the redshift range z w 0.1 — 0.6. The sample is selected 
from images obtained with the Wide Field and Planetary 
Camera 2 (WFPC2) on board the Hubble Space Telescope 
(HST) taken from the HST- Archive. 

The motivations of the project and the first results of 
this study have been discussed discussed in the paper pre- 
senting the analysis of the pilot sample (Treu et al. 1999; 
hereafter T99). In this paper, we describe in detail the se- 
lection of a larger sample of 35 objects, and we present the 
surface photometry and the kinematic spectroscopic mea- 
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surements. An extended discussion of the state of the art 
in this research area and of the objectives of our work is 
given in a companion paper (Treu et al. 2001, submitted to 
MNRAS; hereafter PHI), where the evolution of the Funda- 
mental Plane and Kormendy relation is analyzed. 

This paper is organized as follows. The sample selection 
criteria and the characteristics of the sample are discussed in 
Section ^. In Section ^, we report on the measurement of the 
photometric parameters. The measurement of redshifts and 
central velocity dispersions is described in Section ^| The 
caveats of the photometric and kinematic measurements are 
extensively discussed. In particular, we discuss: i) the accu- 
racy that can be achieved in measuring the combination of 
photometric observables that enters the Fundamental Plane; 
ii) the effects of resolution on the measurement of inter- 
nal velocity dispersion, by means of Montecarlo simulations; 
lii) the stability of the measured internal velocity dispersion 
with respect to the rest frame spectral region used for the 
fit. A summary is given in Section [s| 

The Hubble constant is assumed to be Ho = 
50/i5okm s _1 Mpc _1 . The matter density of the Universe and 
the cosmological constant in dimesionless form are indicated 
as Q an Q\ respectively. 



2 SAMPLE SELECTION 
2.1 Selection criteria 

The targets used in this study have been chosen from a 
sample of random early- type galaxies (E/SO) found in the 
WFPC2 parallel images collected by the Medium Deep Sur- 
vey (Griffiths et al. 1994). The Medium Deep Survey consists 
of a database of 250 WFPC2 fields with at least 1 image in 
each of the filters F606W and F814W, which broadly corre- 
spond to Johnson- Cousins filters V and I; we refer to these 
filters as V606 and Igi4. The images have been reduced and 
analyzed by the MDS group and the complete catalog is 
available to the scientific community (Ratnatunga, Griffiths 
& Ostrander 1999). 

For the sample selection we used the following criteria: 

(i) Morphology clearly denned as early-type. We 

selected the galaxies classified by the MDS as pure bulges or 
disc plus bulge objects with bulge-to-total luminosity (B/T) 
greater than 0.3 (the MDS group considers r 1 / 4 and ex- 
ponential luminosity profiles, indicating them as bulge and 
disc luminosity profile respectively; see the MDS publica- 
tions for details) . All galaxies were inspected by eye to check 
their morphology and some of them were rejected as mis- 
identifications (stars or galaxies with clearly present spiral 
structure). We did not exclude from the sample galaxies with 
nearby companions or disturbed galaxies, in order not to 
bias against merging and interaction. 

(ii) Sufficiently high luminosity: (Isi4 < 19.3), which 
is the effective limit with EFOSC2 at the ESO-3.6m with 
the adopted grism #9 in normally good seeing conditions 
(O'.'8-l"). 

(iii) High galactic latitude: (j&j > 15°), to limit fore- 
ground star contamination and galactic extinction. 

(iv) Low extinction: (E(B-V) < 0.2), as measured by 
Schlegel, Finkbeiner & Davis (1998), to reduce the possibil- 
ity of photometric errors due to patchy reddening. 



(v) Appropriate colour: (0.95 < V606-Isi4 < 1.7), 
which is the entire range spanned by E/S0 in the approxi- 
mate redshift range z = 0.1 — 0.6. 

(vi) Clustering properties. Parallel images are taken 
at random pointings within a few arcminutes of the primary 
target. Therefore, when the primary target is a known clus- 
ter, the parallel image can also contain part of the cluster. 
In order to avoid such contamination, and to obtain an un- 
biased field sample, images centred on known clusters were 
discarded from the analysis. 



2.2 Discussion of the adopted selection criteria 

Apart from the natural requirements on the magnitude limit 
and the limits on galactic latitude and foreground extinc- 
tion, some selection criteria deserve a thorough discussion. 
In fact, understanding the selection criteria is a crucial step 
for a meaningful interpretation of the data (e.g., see the 
discussion in Schade et al. 1999 on morphology and colour 
selection criteria). 

Morphological Selection. We checked our morphologi- 
cal selection criterion (B/T > 0.3) by studying the distribu- 
tion of B/T values for the galaxies morphologically classified 
as early- type by Abraham et al. (1996). By matching abso- 
lute coordinates, we cross-correlated the MDS and Abraham 
et al. catalogs down to \%\a = 21 recovering the B/T value 
assigned by the MDS to each galaxy visually classified as 
an early- type by Abraham et al. (1996). The distribution of 
B /T values found by the MDS group is plotted in Figure |l| 
We conclude that our threshold in B/T does not exclude a 
priori any significant part of the early-type galaxy popula- 
tion. 

Colour selection. Our experimental setup was tuned to 
work in the redshift range z = 0.1 — 0.6. Therefore, we chose 
very loose colour cuts that at the same time could reject 
E/S0 outside our target redshift range (z < 0.1 or z > 0.6) 
and include the widest range of stellar population properties. 
The colour-redshift relation for our sample is shown in Fig- 
ure^. It is noticed that at any given redshift there are objects 
with early-type morphology that are significantly bluer than 
the reddest ones. The effect of the adopted colour selection 
criterion on the study of the evolution of stellar populations 
is discussed in PHI. 

Clustering Properties. We did not attempt to identify 
and to exclude members of groups or poor clusters. In fact, 
four of the six objects for which we have obtained FP param- 
eters in T99 are likely to be members of a group (or a poor 
cluster) . Hence, the early-type galaxies in our sample should 
be representative of a random pointing, magnitude-selected, 
sample of objects. 



2.3 The sample 

We took spectra of 25 galaxies satisfying the adopted colour 
and magnitude selection criteria. For 22 of them the signal- 
to-noise ratio achieved is sufficient to measure accurate ve- 
locity dispersions (see Section^). By rotating the slit on the 
sky, we observed spectra for 10 additional early-type galax- 
ies as secondary targets. For 6 secondary early-type galaxies 
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Figure 1. Comparison between the morphological classification 
of Abraham et al. (1996) and our selection criteria. The catalog of 
all galaxies visually classified as E/SO by Abraham et al. (1996) 
was correlated with the MDS catalog by matching the absolute 
coordinates. The distribution of the bulge-to-total luminosity ra- 
tio (B/T) measured by the MDS for this catalog is plotted. A 
small value of the B/T cut (we adopted B/T> 0.3) is needed in 
order not to reject a significant number of early-type galaxies. 



the signal-to-noise ratio achieved is sufficient to measure ve- 
locity dispersions. The redshift has been measured for all 
the 35 galaxies for which we present photometry. 

In Section ^, where we are mostly interested in dis- 
cussing the measurement of velocity dispersions, we will 
consider all the galaxies including the secondary targets. In 
turn, when analyzing the results of this study in terms of 
evolution of the stellar populations of E/SO (as in PHI), we 
will consider only the galaxies satisfying the adopted selec- 
tion criteria. 

Table [l] lists the observed objects (the digit in the 
adopted name indicates the observing run when the spectra 
were taken, see Section ^) with their coordinates and MDS 
identification. The objects are flagged to indicate whether 
they satisfy all the selection criteria (sel=2), all criteria ex- 
cept that they are fainter than the limit Isi4 = 19.3 (sel=l), 
or they are added only as the best secondary target in the 
field (sel=0). For some of the objects the MDS photometry 
was not available (galaxy D3 is not considered by the MDS 
group because it is close to the edge of the combined frame 
used, see Table ^; sel=nn), or was not reliable (chip edge 
for galaxy Q3, dust lane for galaxy F3; sel=-l). For easy 
reference Table [jj also lists for each object the number of 
spectroscopic observations, the total spectroscopic integra- 
tion times, the S/N achieved, and the redshift of the object. 

Since we wish to use this sample to infer general prop- 
erties of the population of early-type galaxies, it is very 
important to understand and quantify any bias that may 
be introduced in the selection process. The distributions of 
apparent magnitude, colour, and B/T of the galaxies with 
measured velocity dispersion and satisfying the relevant se- 



Figure 2. Observed V606-I814 colour (corrected for foreground 
extinction) as a function of redshift for our sample. The galaxies 
are plotted as pentagons (E), squares (E/SO) or triangles (SO) ac- 
cording to their morphological classification. Morphological classi- 
fication based on HST morphology and surface photometry is dis- 
cussed in Section 3.4. Peculiar galaxies are identified by an open 
symbol. In order of increasing redshift they are: F3 (z = 0.222; 
central dust lane, see Section g), E2, C2 (z = 0.398, 0,425 strong 
emission lines; see the spectrum of E2 in Figure pj), and B3 
(z = 0.490; strong Balmer absorption features, see Figure 
Note that at any given redshift morphologically classified objects 
span a significant range in colour. 



lection criteria are shown in Figure |^, compared to those of 
the galaxies with I 8 i 4 < 19.3 and 0.95 < V606-Isi4 < 1.7 in 
the MDS. As expected from our morphological classification 
we have a higher frequency of objects with large B/T. No 
major bias is found to be present in the colour distribution. 

In addition it is important to notice that, because of 
the sele ctio n effects (Igi4 < 19.3) and resolution limits (see 
Section 



4.3), our sample is limited to relatively massive lumi- 



nous early-type galaxies. The distribution of absolute mag- 
nitudes is plotted in Figure ^ together with the location 
of the characteristic magnitude M* of the Schechter lumi- 
nosity function (Schechter 1976). The bulk of our sample of 
galaxies is comprised between A4* and M.* — 1. 



3 PHOTOMETRY 

The images are taken from the Medium Deep Survey of the 
Hubble Space Telescope (Griffiths et al. 1994). For each 
MDS field there are images from WFPC2 through filters 
V606 and Isi4. Table [| lists total exposure times and num- 
ber of exposures. 

The data reduction procedure has been described in full 
detail in T99 (see also Treu 2001) and the description will 
not be repeated here. 
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Table 1. Galaxy identification. For each galaxy we list the adopted name (gal), the MDS field nam e and identification number within 
the field (MDS ID), MDS coordinates (RA DEC), satisfaction of sel ecti on criteria (sel, see Section [j,3| ), and whether a was measurable or 
not (cr). Galaxy D3 is not present in the MDS catalog (see Section 2.3) and is marked with nn in columns MDS ID and sel. In addition, 
for each galaxy we list: number of spectroscopic observations (n), total spectroscopic exposures time (texp in seconds), and average S/N 
per pixel (S/N). Finally, the redshift is listed in the last column, with four digits if determined with the GHFF (see Section 4), with 
three if determined by recognizing the most important absorption features. 

gal MDS ID RA(J2000) DEC(J2000) sel a n texp S/N z 



Al 


ut800.1 


14:42:14.72 


-17:11:43.7 


2 


y 


4 


7200 


65 


0.1466 


Bl 


ut800.8 


14:42:12.87 


-17:11:57.8 





y 


4 


7200 


40 


0.146 


Dl 


ur610.6 


13:25:40.81 


-29:54:13.5 


2 


y 


4 


12000 


22 


0.3848 


El 


U5405.62 


17:57:34.29 


+04:51:32.8 


2 


y 


5 


13200 


17 


0.2940 


Fl 


U5405.91 


17:57:34.58 


+04:51:42.4 


2 


n 


5 


13200 


12 


0.2935 


Gl 


U5405.18 


17:57:33.54 


+04:51:42.1 


2 


y 


3 


10800 


35 


0.2946 


11 


U5405.47 


17:57:33.47 


+04:52:07.3 


2 


y 


3 


10800 


22 


0.2929 


01 


ustOO.15 


10:05:47.23 


-07:42:05.9 





ii 


10 


36000 


7 


0.647 


PI 


ust00_16 


10:05:47.89 


-07:41:53.1 


1 


ii 


10 


36000 


6 


0.550 


B2 


upuOO.13 


09:59:18.54 


-22:54:46.5 


1 


ii 


2 


7200 


6 


0.463 


C2 


upzOCLlO 


13:22:09.46 


-36:40:10.9 


2 


ii 


2 


7200 


9 


0.425 


D2 


urwl0_3 


11:21:24.02 


-24:55:18.6 


2 


y 


2 


7200 


33 


0.1923 


E2 


uu301.1 


15:56:33.17 


+11:08:25.9 


2 


y 


i 


13500 


21 


0.3977 


F2 


uu301.5 


15:56:32.86 


+11:07:57.6 


2 


y 


4 


13500 


18 


0.3364 


A3 


uxlOO.2 


14:45:09.52 


+10:01:37.5 


2 


y 


3 


7200 


19 


0.2311 


B3 


uxlOO.5 


14:45:09.76 


+10:02:09.7 


2 


n 


3 


7200 


12 


0.490 


C3 


ud200.2 


21:58:40.34 


-30:22:32.4 


2 


y 


2 


3600 


27 


0.1057 


D3 


uha01_nn 


21:32:32.79 


+00:15:04.8 


nn 


y 


2 


3600 


16 


0.1511 


E3 


uha01_12 


21:32:34.28 


+00:14:13.8 


2 


y 


4 


10800 


17 


0.2631 


F3 


uhaOl-10 


21:32:33.89 


+00:14:36.1 


-1 


n 


2 


7200 


10 


0.2220 


G3 


ubz09_2 


00:49:35.21 


-52:03:58.7 


2 


y 


7 


23400 


21 


0.4081 


13 


ucilO.2 


01:24:39.56 


+03:52:21.4 


2 


y 


7 


25200 


22 


0.4103 


L3 


ucilO.8 


01:24:38.07 


+03:52:06.9 


1 


y 


7 


25200 


13 


0.4141 


M3 


uvlOO.8 


15:06:20.65 


+01:44:02.6 


2 


y 


4 


7200 


16 


0.3374 


Q3 


uvlOO.l 


15:06:21.45 


+01:45:16.7 


-1 


y 


2 


3600 


11 


0.1139 


R3 


ua-00.2 


01:02:24.52 


-27:10:39.9 


2 


y 


4 


12400 


24 


0.3482 


S3 


ua-00.3 


01:02:22.78 


-27:10:58.7 


2 


y 


5 


16000 


18 


0.3559 


T3 


ucsOl.2 


02:56:20.69 


-33:21:22.3 


2 


y 


2 


1800 


14 


0.1174 


U3 


ucs01_4 


02:56:25.11 


-33:22:08.1 





y 


4 


7200 


19 


0.2401 


A4 


ucs01_3 


02:56:25.26 


-33:21:24.7 


2 


y 


2 


3600 


25 


0.1171 


131 


ua400_3 


00:24:56.52 


-27:16:10.2 


2 


y 


2 


5356 


19 


0.2112 


D4 


ufrOO.2 


04:08:38.28 


-24:23:53.3 


2 


y 


5 


13567 


32 


0.2848 


E4 


ufr00.3 


04:08:38.94 


-24:25:44.4 





y 


5 


13567 


27 


0.1851 


H4 


uim03_l 


03:55:33.99 


+09:43:11.5 


2 


y 


7 


25200 


25 


0.3399 


14 


uim03_7 


03:55:35.09 


+09:43:32.6 


2 


y 


7 


25200 


19 


0.3369 



For each galaxy we fitted isophotal profilesQ and 2D 
models convolved with the Tiny Tim PSF (Krist 1994; see 
T99). We define the photometric parameters, effective radius 
re, and effective surface brightness SBe, as the best-fitting 
parameters of the r 1 / 4 component. As in T99, we will denote 
the angular effective radius (in arcseconds) by re and the 
linear effective radius (in kiloparsecs) by Re- 

The isophotal profiles were fitted with an r 1//4 profile, 
an r 1,/4 plus exponential profile, and an exponential pro- 
file. All galaxies in the sample were best-fitted by the r 1//4 
or r 1//4 plus exponential luminosity profile. Fifteen galaxies 
were best-fitted by a pure r 1//4 profile. Fifteen gave insignif- 
icantly smaller residual when an exponential component was 
added (and negligible changes in the r 1//4 photometric pa- 
rameters). For five galaxies the addition of the exponential 
component improved significantly the fit (see Section 3.4) 
and changed the parameters. In the following we adopt the 



* Surface brightness as a function of circularized radius, as de- 
fined in T99. 



values obtained by fitting a pure r 1//4 profile, in all cases ex- 
cept for the five galaxies for which the exponential is needed 
(galaxies Bl, D2, E3, S3, A4), where the two components 
are fitted and the parameters of the r 1//4 component are 
adopted. The 2D models were pure r 1 / 4 luminosity pro- 
files, and the fit was performed with two weighting schemes, 
least squares (2Dls) and least \ 2 (SDhv 2 ). 

The larger sample available for the study presented in 
this paper, with respect to T99, allows us to perform a sta- 
tistical analysis of the photometric results. In the following 
we will give relative errors in per cent when referred to ef- 
fective radii, in magnitudes when referred to magnitudes or 
surface brightnesses, and in absolute terms when referred to 
logarithmic quantities. 

The values obtained from luminosity profile fitting and 
2D1^ 2 have an average difference of 15 per cent in re (in 
F606W; 14 per cent in F814W) and r.m.s. scatter of 16 per 
cent (F606W and F814W), if we restrict the study to the 
objects without a significant disky component. If the five 
objects with discs are included the r.m.s increases to 53 per 
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Figure 3. Comparison of the properties of the galaxies in our 
sample with measured velocity dispersion to the ones of the galax- 
ies in the full MDS satisfying our colour and magnitude selection 
criteria (Igi4 < 19.3, 0.95 < V606-I814 < 1-7). From top to bot- 
tom, we plot the distribution of bulge-to-total luminosity ratio 
(B/T), magnitude (1814) and colour (V6O6-I814) f° r °ur sample 
(points with error bars, computed using Poisson statistics) and 
for the MDS catalog before and after B/T cut (thin and thick 
histograms). We note that, as expected since we avoided early- 
spirals, we have a marginally higher frequency of large B/T values 
with respect to the entire catalog, and the small B/T values are 
rejected. Furthermore, we have a smaller frequency of objects in 
the faint magnitude bin, while no significant bias is present in the 
colour distribution. 



Figure 4. Distribution of absolute magnitude for the galaxies 
satisfying the colour and magnitude selection criteria (empty his- 
togram). The subset with measurable a (S/N > 12) is plotted as 
hatched his togram. The subset with a above resolution limits (see 
Section |4.3| ) is plotted as filled histogram. The distributions for B 
and V rest frame luminosities are shown in the upper and lower 
panel. The thick vertical dashed line represents the characteristic 
magnitude M* as measured by Marzke et al. (1998) in the local 
Universe (in the V band M* is obtained by assuming B-V = 0.95, 
the colour of a single-burst stellar population, age of 12 Gyr, so- 
lar metallicity, and Salpeter IMF (Salpeter 1955), computed using 
Bruzual & Chariot 1993, GISSEL96 version). The values of the 
cosmological parameters Q = 1 and Qa = are assumed. 



cent (F606W) and 37 per cent (F814W); we emphasize that 
the 2D models have pure r 1 ^ 4 luminosity profiles and there- 
fore an increased scatter is expected, because of the different 
modelling. The least square fitting technique produces larger 
differences and scatter in the results. Compared to the pro- 
file fitting method, the average difference in re is 30 per cent 
(F606W) and 28 per cent (F814W) with an r.m.s. scatter of 
20 per cent (F606W) and 18 per cent (F814W). The r.m.s. 
scatter increases to 32 per cent and 27 per cent if the galaxies 
with a significant disc component are considered. 

As noted by many authors, the effective radius and ef- 
fective surface brightness are correlated observables and the 
combination that enters the FP, 



FP ph =\ogr e - PSB e , 



(1) 



is particularly robust (see the careful discussion by Kelson 
et al. 2000 of the results obtained by fitting a variety of mod- 
els, including the Sersic r 1//n profile, with a variety of tech- 
niques, including growth curves, luminosity profile fit and 
2D fit). The values of the quantity^ FP p h obtained with the 
luminosity profile fit and 2Dlx 2 differ on average by -0.003 



(F606W) and 0.001 (F814W) with an r.m.s. scatter of 0.037 
(F606W) and 0.034 (F814W), by including all the galax- 
ies. In addition, FP p h is also remarkably stable for 2Dls; 
the average difference with respect to the results obtained 
by fitting luminosity profiles is -0.029 (F606W) and -0.027 
(F814W) with an r.m.s. scatter of 0.031 (F606W) and 0.026 
(F606W). 

Because of the higher scatter and difference obtained 
with the 2Dls method, we adopt as best estimate of the 
photometric observables the average of the results obtained 
with the luminosity profile and the 2Dlx 2 method. The best 
estimates of the photometric observables are listed in Ta- 
ble together with the uncertainty taking into account er- 
rors related to sky subtraction, flatfielding and fitting tech- 
nique, computed as described in Section 3.3. As a check, we 



compared the magnitudes we derived here to those derived 
by the MDS group: the average difference and r.m.s scat- 
ter are 0.02, 0.15 (V 6 oe), 0.01, 0.12 (Ism), and 0.01, 0.11 
(V606— 1814), after excluding galaxy F3 (see below). The 
comparison takes into account the different photometric zero 
point s adopted by the MDS group (s ee the MDS web-site at 
URL http://archive.stsci.edu/mds/) and the ones given on 



t The values listed here are computed using j3 = 0.32 (F606W) 
and j3 = 0.328 (F814W). However, since the slope j3 is very well 
determined (see Pahre, Djorgovski & de Carvalho 1998) and varies 



very little with wavelength, the estimate of the errors changes 
negligibly within the observed range of values for /3. 
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Table 3. Observed photometric structural parameters. The uncertainty on the combination that enters the FP, FP V ^ = log Re — PSBc, 
is listed as SFP. The values of SFP listed here are calculated using /3 = 0.32 (F606W) and /3 = 0.328 (F814W; see text). The uncertainty 
on the colour is listed as 568. 
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eo 




CO 




<5FP 6 




1-814 




Co 






eo 
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Al 


17 


.35 ±0.06 


20 


.39 ±0.06 


1 


.62 ±0. 


.08 


0.011 


16. 
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the HST User Handbook that are adopted in the present 
paper. 

Galaxy F3 required a special treatment, which is worth 
an extended discussion. The central part of the galaxy hosts 
a dust lane, which alters significantly the luminosity profile. 
In order to recover the underlying surface luminosity we 
performed the following correction. Assuming that the dust 
lane effects can be neglected at large radii and that colour 
gradients are negligible for the correction, we computed the 
observer frame extinction map, 

E6S(x, y) = ix 6 (x, y) - us {x, y) - M68 (oo) (2) 

where ns(x, y) and (i&(x, y) are the surface brightness at any 
given pixel through filters F606W and F814W and /i6s(oo) is 
the colour at large radii. Based on the extinction law given 
by Cardelli, Clayton & Mathis (1989) and on the redshift 
z = 0.222, we calculated the relations between the extinction 
map and the absorption in the individual filters (Aq and Ag): 

A 6 {x, y) = 3.401_E68(:r, y) (3) 

As(x,y) = 2A01E68(x,y). (4) 

The results of the correction are good as it can be judged 



from the quality of the luminosity profile shown in Figure |5| 
and of the images shown in Appendix [X]. However, because 
of this correction, the galaxy is flagged by an open symbol 
in Figure bl 



3.1 Galactic extinction 

We used E(B-V) values given by Schlegel et al. (1998) and 
the relations A 6 = 2.889 E(B-V) and A s = 1.948 E(B-V) 
calculated by Schlegel et al. (1998). The galactic extinction 
for each field is listed in Table ^. We compared the extinc- 
tion given by Schlegel et al. (1998) with the values given by 
Burstein & Heiles (1982). The average difference in E(B-V) 
is 0.017, consistent with the different zero point used by the 
authors (see Schlegel et al. 1998), and the scatter is 0.012. 



3.2 Rest frame photometric quantities 

As in T99, the rest frame photometric quantities are com- 
puted from the observed ones using the K-color correction 
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Figure 5. Luminosity profiles, I. The luminosity profiles obtained by isophotal fitting are plotted together with the best-fitting r 1 / 4 model 
convolved with the PSF. Surface brightness (fi) is plotted in mag arcsec -2 , radius (r) in arcsec. The upper curves represent luminosity 
profiles through filter F8f4W, the lower through filter F606W. The profiles of galaxies Al, Bl, Dl, II, Gl, El, Fl, Ol, PI are taken 
from T99. Galaxy Bl shows a spiral pattern in the residuals of the two-dimensional fit (see T99). Galaxy D2 is fitted by an r 1//4 plus 
exponential model. 



(Amvg, Amse; listed in Table |^) defined as: 



Arrive 



-2.5 log 



/ F x (X)S v (X)dX 



+2.5 log 



/ F X (X)S 8 (X(1 + z))dX 



Z P 8 > ( 5 ) 



Am B 6= -2.5 log J F x (X)S B (X)dX + zp 



+2.51og jF x (X)S 6 (X(l + z))d\ -zp 6 , (6) 



where zp and S(X) are the zero points and transmissions of 
the band passes in the Landolt system, while F\(X) is the 
flux per unit wavelength of the model spectrum. The model 
spectrum used is the synthetic spectrum (from Bruzual & 
Chariot 1993; GISSEL96 version) of a single-burst stellar 
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Figure 6. Luminosity profiles, II. As in Figure H. Galaxy F3 hosts a dust lane in the ce ntra l region; the profile shown here is corrected 



for extinction as described in Section ^. Galaxy E3 has a clear spiral pattern, see Figure 



Al 



population that bes t reproduces the observed colour (see 
T99 and Section E|). 



The radii r g g and r g y are calculated for the central 
rest wavelength of each filter B and V by linear interpo- 
lation/extrapolation in wavelength between the measured 
radii in the WFPC2 niters F606W and F814W at their ob- 
served central wavelengths. We take the central wavelengths 
of B, V, F606W, and F814W to be 4400, 5500, 5935, and 
7921 A, respectively, resulting in the following formulae for 



the angular sizes: 

ree = {re F6 06w[7921 - (1 + 2)4400] 

+r eF 8i4w[(l + z)4400 - 5935]}/1986 (7) 
rev = {reF606w[7921 - (1 + z)5500] 

+r eF 8i4w[(l + 2)5500 - 5935]}/1986 . (8) 
The rest frame surface brightness is computed as: 
SBev = SBe8-101og(l + z)+Am V8 -A 8 (9) 
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Figure 7. Luminosity profiles, III. As in Figure]^ Galaxies S3 and A4 are fitted by an r^/^plus exponential model. 



SB es = SB e6 - 101og(l + z) + Am B6 - A 6 



(10) 



Table q lists for each galaxy the effective radius in kpc, com- 
puted for a cosmological model with Q = 1, = 0, hso = 1, 
and the effective surface brightness. 



3.3 Error analysis 

The main sources of error are given by: 

(i) Fitting technique. We estimate the error in terms of 
the half-difference of the values found with the two indepen- 



dent fitting techniques used (luminosity profiles and 2D1% 2 
fitting). 

(ii) PSF modelling. We estimate this effect by fitting a 
small galaxy in the corner of the chip (D3) and a large one 
in the centre of the chip (H4) with a set of different PSFs. 
We use a Tiny Tim PSF generated on the galaxy position, 
a star present on the chip, and a Tiny Tim PSF generated 
on the position of the star. For D3, which is observed on 
three frames with different subpixel pointings, we produce 
a PSF by summing three subsampled PSFs centred on the 
corresponding position within the pixel. The r.m.s. scatter 
of the effective radius derived with the various PSFs is 3 per 
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Table 5. Rest frame photometric structural parameters. Effective radii are in hT^ kv c (for Q = 1 and f2yv = 0). The error on the rest 
frame colour B-V is listed as <5BV. The morphological classification (Class; Section |3.4| ) is a lso listed. Photometric structural parameters 
for SO and S refer to the r 1 / 4 component. Errors are computed as described in Section p.3|. 
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cent for D3 and 2 per cent for H4 and the scatter of FP p h is 
respectively 0.007 and 0.004. Therefore, this is a negligible 
source of error for our purposes. 

(iii) Sky subtraction and flat fielding. We measured 
the contribution to the error given by these uncertainties 
by varying the background by 1 per cent (which we take as 
upper limit to the errors in sky subtraction and flat fielding) 
and by running the 2D fits. These errors are smaller than 
the ones due to the fitting technique, but in some cases they 
are not negligible. 

(iv) Galactic extinction. We adopt the error estimate 
by Schlegel et al. (1998) for the reddening correction <5E(B- 
V)=0.16E(B-V). 

(v) K-color correction. For each galaxy, we com- 
puted the colours of single-burst stellar population spectra 
(Bruzual & Chariot 1993; GISSEL96 version) with ages be- 
tween 0.5 and 20 Gyr, redshifted to the galaxy redshift. We 
also computed K-color corrections for each synthetic spec- 
trum. By interpolation we found the K-color correction cor- 
responding to the colour of each galaxy. This procedure ap- 
plied to the colours shifted by a standard deviation provided 
upper and lower limits to the K-color correction. The error 
on the K-color correction is the half-difference of the up- 



per and lower limits. Table g lists the K-color corrections 
(Am_B6 and Amyg) together with their errors (5Ams6 and 
SA.mva)- 

The total error on the measured parameters listed in Table ^ 
is the quadratic sum of terms i and iii. The uncertainties 
on the extinction correction and on the K-color correction 
are added in quadrature to produce the total errors listed 
in Table ^. The systematic uncertainty in the zero point 
is not considered here and will be taken into account when 
comparing the intermediate redshift results to the local ones 
in PHI. 



3.4 Classification 

The objects are classified in a simplified scheme (E, E/S0, 
SO, S) based on their morphology and luminosity profile. 
Among the early-type systems, objects that are best-fitted 
by a pure r 1 ^ 4 profile are classified as E; objects for which 
the addition of an exponential component gives a slightly 
better fit with no significant change in the structural pa- 
rameters are classified as E/S0; objects that are significantly 
better described by an r 1 ^ 4 plus exponential profile are clas- 
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Table 2. Photometric data: number of exposures (in each of the 
two filters, n6 and ns) and total exposure times (in each of the 
two filters, texp6 and texps, in seconds). The fields uhaOO and 
uhaOf and uimOl and uim03 overlap partially. Some of the fields 
have been imaged with large dithering and therefore the exposure 
time is not uniform across the field. When two galaxies in the same 
field have been imaged with different exposure times we list the 
field twice, specifying the two exposure times. In column note we 
specify which galaxy has been imaged with the listed exposure 
times. Column E(B-V) lists the galactic extinction as given by 
Schlegel et al. (1998). 
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sified as SO. Objects that show a spiral pattern (Bl, E3) are 
classified as S. The classification is listed in Table ^| Note 
that in this classification the usage of the symbol E/S0 is 
different with the respect to the rest of the paper, where it 
indicates the entire class of early- type galaxies. 



4 SPECTROSCOPY 

A total of 13 dark nights at the ESO 3.6m telescope were 
awarded for this spectroscopic project. 

The first pilot observing run (run 1, described in T99) 
took place in 1996, April 18-22. Long slit (l'/5 slit width) 
spectroscopy was obtained with the EFOSC spectrograph, 
with the Orange 150 grism in the spectral range between 
approximately 5200 and 7000 A (512 x 512 CCD, with pixel 
size 30/im equivalent to about 3.5A x 0'.'57). The estimated 
seeing, measured with R-band acquisition images, was 0'.'8 
FWHM. 

The following runs (run 2, 3, 4) took place in 1999, 
March 19-21, August 11-14, and November 2-6. The imager- 
spectrograph EFOSC2 was used with grism # 9, which cov- 
ers the range 4700-6700 A, with pixel size of approximately 
2AxO'.'32 in 2 x 2 binning mode. A narrower slit (1") was 
used in order to achieve a better resolution. The seeing 
ranged from 0'.'8 — '.'2 during run 2 and 3. Run 4 was subject 
to highly variable weather conditions: two nights were lost 
for clouds and bad seeing (2 — 3"), one night had average 



Tab le 4. K-color correction, computed as described in Sec- 
tion E rror s ((5Amgg and <5Amyg) are computed as described 
in Section p. A. 
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seeing conditions (0'.'8 — 1"), and during the last night we 
had some episodes of excellent seeing (image quality 0'.'6). 

He-Ar lamps in the pointing position were taken be- 
fore and after every set of 2 exposures to provide accurate 
wavelength calibration. Standard stars were observed at the 
beginning and at the end of every night to provide relative 
flux calibration. 

The spectra were reduced as described in T99 (see also 
Treu 2001). For every galaxy, Table |l| gives total integration 
time, number of exposures, S/N achieved. In Figures ^ to ^ 
the spectra with highest signal-to-noise collected are shown. 
The run is identified by the digit in the name. The spectra 
taken during run 1, shown in T99, are not repeated here. 

For runs 2, 3, and 4, by taking advantage of the higher 
number of spectra available with the same instrumental 
setup, we were able to determine the resolution of our final 
combined spectra more accurately than in run 1 (see T99). 
We selected a set of 3 sky lines and 12 He-Ar lines that are 
unblended at our resolution and sufficiently strong and iso- 
lated and we measured their FWHM on each combined sky 
and He-Ar calibration spectrum. The results were very sta- 
ble from frame to frame and from sky to He-Ar lines, so we 
fitted the entire set of measurements as a function of wave- 
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Figure 8. Wavelength calibrated spectra of the galaxies with average S/N per pixel greater than 10 taken with EFOSC2, I. The shaded 
spectral regions were affected by strong sky-emission lines. These regions have not been used in the kinematic fit. The standard deviation 
per pixel is also shown (lower curves). The Calcium doublet (Ca H 3968.5 A and K 3933.7 A), the G band (4304.4 A), H/3, and Mg,, 
(5175.4 A) are the most prominent absorption features in this spectral region. Note the extremely strong Balmer absorption features in 
the spectrum of B3 and the NaD (5892.5 A) line in the spectrum of the lowest redshift objects such as C3. Note [OIIJ3727 (A)in emission 
in the spectrum of E2. The central velocity dispersion in kms -1 is shown in the upper left corner when measured (see Section 4). 



length. The measurements of the resolution (a s in km s _1 ) 
are well fitted by a parabola, as shown in Figure [Hj, left 
panel. 



4.1 Fitting kinematic parameters 



The internal velocity dispersion is derived by comparison of 
the observed galaxy spectrum to a stellar spectrum taken 
with the same resolution under the assumption that the 
galaxy spectrum can be well described by the sum of the 
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Figure 9. Wavelength calibrated spectra of the galaxies with average S/N per pixel greater than 10 taken with EFOSC2, II. See caption 
to Figure ^. 



spectra emitted by stars with a Gaussian velocity distribu- 
tion. 

In the first exploratory study (T99) we derived kine- 
matic information for our sample of galaxies with two dif- 
ferent and independent methods: i) the Gauss-Hermite 
Fourier Fitting Software (Franx, Illingworth & Heck- 
man 1989; van der Marel & Franx 1993), hereafter GHFF, 
and ii) the Fourier Quotient (Sargent et al. 1977) in 
the version modified by Rose (see Dressier 1979) and by 
Stiavelli, M0ller & Zeilinger (1993). The results of the two 
methods were mutually consistent, within the estimated er- 



rors. However, the GHFF allows for a reliable estimate of 
the errors, is less sensitive to template mismatches (van der 
Marel & Franx 1993), and provides insight in the fitting pro- 
cedure. Finally, it provides a quality parameter (the \ 2 ) and 
the residual spectrum. Therefore, here we will use only the 
GHFF software. 



4.2 Stellar templates 

A good library of stellar spectra is a key ingredient to mea- 
sure the internal kinematics of a galaxy. It has to cover the 
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Table 6 ■ S tars observed as stellar templates (SK library; see 
Section 4.2 for details). For each star we list the spectral type, 
name of the star, magnitude and colour, taken. 



5000 5500 6000 6500 

A (A) 
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Figure 10. Left: resolution (<r s ) as a function of wavelength as 
determined from the width of sky and He-Ar lines for the EFOSC2 
setup. Right: resolution (<r t ) as a function of wavelength for the 
SK stellar template library (the points represent the measured 
resolution, t he s olid line is the best fitting parabola; see Tabic ^ 
and Sec tion |4.2| ); the resolution of the JHC library (4.76A see 
Section 4.2) is shown for comparison as a dashed line. 



largest possible range of spectral types so that stellar pop- 
ulation effects can be studied. In addition, the spectral res- 
olution of the stellar template (at) must be better than the 
instrumental resolution of the setup used for galactic spec- 
troscopy. The stellar library used in T99 (Jacoby, Hunter & 
Christian 1984; hereafter the JHC library), given its resolu- 
tion (~ 4.5 A FWHM, i.e. a t = 133 km s _1 at 4300 A), can 
only be used out to z ~ 0.2 with our spectroscopic setup. 

David Soderblom and Jeremy King kindly provided us 
with a library of stellar spectra at higher resolution for the 
highest redshift part of the sample (hereafter the SK li- 
brary). The stars, listed in Table [| were observed at the 
Kitt Peak National Observatory Coude Feed Telescope, with 
the Coude CCD Spectrograph, on May 22 1999. The stars 
were chosen according to observability from a list of giant 
stars covering all spectral types from AO to K7. In addition, 
a few dwarf star spectra were also taken. The spectral range 
covered is 3120-5510 A. The resolution of the SK library, 
as measured from calibration lamps, is shown in Figure hoi 
right panel. 

Unfortunately, the limited red coverage of the spectra 
in the SK library did not allow us to use them for the low 
redshift objects of the sample (z < 0.21). Therefore, we 
used the JHC library for the low redshift part of the sample 
(z < 0.21) and the SK library for the high redshift part of 
the sample (z > 0.21). 

In order to derive an accurate and consistent value for 
the resolution of the JHC library, we measured it using the 
SK library. We divided the overlapping region (3510-5510 A) 
into four intervals of 500 A each and we fitted the "velocity 
dispersion" of each interval of the JHC spectra using the 
SK template that more closely matches the spectral type. 
By adding in quadrature the resolution of the SK library, we 
found the intrinsic resolution of the JHC library (4.76±0.30 
A FWHM), in agreement with the value given by Jacoby et 
al. (1984). We used this value as intrinsic resolution of the 
JHC library. 

The stellar templates were redshifted to the redshift of 
each object and then broadened with a Gaussian character- 
ized by a width equal to the quadratic difference between 
the instrumental resolution (a s ) and the resolution of the 
template library (at). 
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4.3 Results and error discussion 

As in T99, we fitted each galaxy with each of the stellar tem- 
plates available^. This procedure allows us to estimate the 
uncertainty given by template mismatches. As noted in T99, 
the late G to early K giant stars provided the best x 2 > the 
smallest residual, and the most stable results. Cooler stars 
produced systematically higher velocity dispersions and hot- 
ter stars systematically lower ones. 

As best estimate of the velocity dispersion and redshift 
(°ghff ' z ghff m Table , we adopted the average of the 
values obtained with templates in the spectral range G4- 
K0III. In order to achieve homogeneous results with the 
two libraries, we used the same number of templates and 
closely matching spectral type distributions: 2xG4, lxG5, 
2xG8, lxG9, lxKO for SK and lxG5, 2xG6, lxG7, lxG8, 
lxG9, lxKO for JHC. In Table ^ we also list the best-fitting 
template for each galaxy. Note that the best-fitting template 
spans the entire range of spectral types. Two estimates of 
the errors are also given, as in T99: the random compo- 
nent, which is the fit formal uncertainty obtained for the 
best-fitting template (S); the systematic component (A), re- 
sulting from template mismatches, computed as the r.m.s. 
scatter of the values obtained with different templates. 

It is intuitive and well known (e.g. Bender 1990; 
J0rgensen, Franx & Kjaergaard 1995; T99) that the S/N 
needed to measure velocity dispersions properly scales with 
the ratio between the instrumental resolution and the veloc- 
ity dispersion itself. In other words, when the velocity dis- 
persion is smaller than the instrumental resolution a higher 
minimum signal-to-noise is needed with respect to the case 



t The spectra from T99 were fitted again using the SK library, 
except for galaxy Al that was fitted with the JHC library spectra 
at the new slightly different resolution. 
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Table 7. Kinematic results from the GHFF fit. For each galaxy we list the average velocity dispersion (fghfT ) together with the random 
error (<5o"g n jj ) and the systematic error (A<Tg n g). The random error is the fit formal uncertainty for the best-fitting template, while 
AtTgjjjf is the scatter of the values obtained with the templates in the range G4-K0 III (see Section 4.3 ). We also list the redshift (.Zghfr ) 
with its errors, the spectral type of t he best- fitting template (Type), the signal-to- noise ratio (S/N), and the matching to the relevant 
selection criteria (sel; see Section 2.3). Only galaxies with S/N > 12 are listed (see Sect ion J4.3| ). The values of velocity dispersion found 
for Bl, D2, E3, Q3, B4, E4 arc too small to be acceptable (see discussion in Section |4,3| ): these objects have been excluded from the 
table. 



gal 


°"ghff 


5a ghS 


Af7 ghff 


2 ghff 


5 z ghff 


A z ghS 


Type 


S/N 


sel 


Al 


227 


16 


14 


0.1466 


4.3-10~ 5 


1.2-10" 4 


G8 


65 


2 


Dl 


290 


21 


12 


0.3848 


6.8-10~ 5 


1.4-10" 4 


G4 


22 


2 


El 


198 


25 


7 


0.2940 


8.4-10" 5 


1.3-10~ 4 


K0 


17 


2 


Gl 


220 


17 


10 


0.2946 


5.9-10" 5 


1.3-10" 4 


G8 


35 


2 


11 


190 


20 


4 


0.2929 


6.5-10~ 5 


1.3-10" 4 


Go 


22 


2 


E2 


219 


12 


3 


0.3977 


4.4-10" 5 


1.4-10~ 4 


G4 


24 


2 


F2 


237 


11 


14 


0.3364 


5.3-10~ 5 


1.3-10~ 4 


G9 


18 


2 


A3 


189 


15 


4 


0.2311 


4.9-10 -5 


1.2-10 -4 


Go 


19 


2 


C3 


146 


15 


9 


0.1057 


3.8-10~ 5 


1.4-10~ 4 


G8 


27 


2 


D3 


136 


18 


8 


0.1511 


4.8-10~ 5 


1.0-10~ 4 


GO 


17 





G3 


244 


11 


7 


0.4081 


4.5-10" 5 


1.4-10~ 4 


Go 


21 


2 


13 


194 


11 


8 


0.4103 


4.M0~ 5 


1.4-10" 4 


Go 


22 


2 


L3 


193 


18 


15 


0.4141 


6.8-10" 5 


1.4-10" 4 


Go 


13 


1 


M3 


228 


15 


5 


0.3374 


5.4-10" 5 


1.4-10" 4 


G4 


16 


2 


R3 


228 


11 


5 


0.3482 


4.0-10~ 5 


1.3-10~ 4 


Go 


24 


2 


S3 


216 


19 


5 


0.3559 


7.2-10~ 5 


1.M0 -4 


Gl 


18 


2 


T3 


133 


23 


16 


0.1174 


5.2-10~ 5 


1.0-10~ 4 


G6 


14 


2 


U3 


238 


14 


21 


0.2401 


5.2-10~ 5 


1.2-10" 4 


GO 


18 





A4 


167 


15 


8 


0.1171 


4.M0~ 5 


1.2-10~ 4 


G7 


25 


2 


D4 


224 


15 


11 


0.2848 


5.1-10" 5 


1.2-10" 4 


Gl 


32 


2 


HI 


217 


12 


5 


0.3399 


4.4-10" 5 


1.4-10" 4 


GO 


25 


2 


14 


243 


20 


13 


0.3369 


7.2-10~ 5 


1.3-10" 4 


G4 


19 


2 



when the velocity dispersion is larger than the instrumental 
resolution. 

In practice, for any given S/N and resolution, there is 
a lower limit on the velocity dispersion measurable without 
introducing significant bias (for example, Bender 1990 con- 
siders half of the instrumental resolution to be the minimum 
measurable value). Quantitatively, these numbers depend on 
the instrumental setup (resolution, sampling). Therefore, we 
addressed the problem by numerical simulation of our mea- 
surement procedure. For a range of values of velocity disper- 
sions (cr in = 65, 75, 100, 125, 150, 175, 200 km s" 1 ) and of 
S/N (10, 12, 15, 18, 25), we created 1000 toy galaxy spectra 
with artificial noise and our instrumental sampling and res- 
olution. We then recovered the velocity dispersion with the 
GHFF software. The results are plotted in Figure [ll]. 

For each input value of <r m we plot the average recov- 
ered o as a function of S/N, using the same stellar spectrum 
as "galaxy" and template (G8 III; panel a) or using a differ- 
ent star as template (G9 III; panel c). It is noticed that the 
smallest values of a m are systematically underestimated by 
the fit for low values of S/N (by as much as 13 per cent for 
o- in = 65 km s _1 and S/N = 10) but for cr m > 100 km s~ x 
the effect becomes negligible (< 2 per cent at S/N = 10, 
< 1 per cent at S/N = 12 for <r m = 100 km s _1 and us- 
ing a different template) . If - instead of the average of a — 
the average of the logarithm of a weighted on the formal 
variance is used as estimator (j0rgensen et al. 1995), the ef- 
fect is reversed: the velocity dispersion is now overestimated 
for small values of <r m (by up to 15 per cent at <r m = 65 
km s _1 and S/N = 10). In addition, high values of a m (from 
200 km s _1 in our simulation) become underestimated by 



2 per cent at S/N = 10. The conclusion is that the system- 
atic uncertainty on small values of S/N is hard to correct, 
so that it is suggested that these data should be rejected. 
In our particular case we will only accept for measurement 
spectra with S/N > 12 (thus rejecting Fl, B3, F3). In Fig- 
ure [ll] we also plot the relative scatter around the mean for 
the same configurations (template G8 III, panel b; template 
G9 III, panel d). For a given S/N the uncertainty on the 
velocity dispersion increases dramatically for the smallest 
values of <r m , while it tends to become constant for velocity 
dispersions much larger than the resolution. This is consis- 
tent with what we find for galaxies B4 and E3 (cghfj = 64 
and 98 km s _1 ): despite the good S/N of the spectrum the 
formal uncertainty on a m is ~ 35 and ~ 18 per cent. We do 
not consider these values to be successfully measured. Simi- 
larly, the small value of velocity dispersion of galaxy Bl was 
rejected in T99. 

Since most of the objects at z < 0.21 have high S/N 
spectra, the problems discussed above are negligible. Nev- 
ertheless, the uncertainty in the JHC spectral library reso- 
lution (~ 120 ± 8 km s _1 at 5000 A) can cause systematic 
effects on the low velocity dispersion objects. We ran ex- 
tensive simulations (as above) to measure this effect. For 
velocity dispersions smaller than the resolution of the li- 
brary, the uncertainty is important (for example ~ 10 per 
cent at 100 km s" 1 with S/N = 25), but it drops quickly 
with increasing velocity dispersion (6 per cent at 125 km s _1 
with S/N = 25; 4 per cent at 150 km s _1 ). In practice, at 
a m = 125 km s _1 this effect is already smaller than the typ- 
ical uncertainty. We will not consider galaxies D2, Q3, E4 
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Figure 11. Results of Montecarlo simulations for the EFOSC2 
setup. For a range of values of velocity dispersion (o" m ) the ve- 
locity dispersion recovered is plotted as a function of S/N (using 
the same star as "galaxy" and template in panel a, and using a 
different star as template in panel c). The scatter of the distribu- 
tion is plotted in panels b (same star as template and galaxy) and 
d (different star as template and galaxy). At S/N < 15 velocity 
dispersions below 100 km s _1 are systematically underestimated. 
Moreover, for small velocity dispersions, the uncertainty is large 
even for higher values of S/N. 



as successfully measured because their velocity dispersion is 
smaller than this limit. 

The uncertainty in the measurement of the instrumen- 
tal resolution of EFOSC (AA = 8.39 ± 0.26A) introduces a 
systematic error in the measurement of the velocity disper- 
sion that is negligible in our case. For example, let us con- 
sider Al, the galaxy with the highest signal-to-noise spec- 
trum. The random error is estimated to be 6 per cent, the 
systematic error due to template mismatches is 7 per cent; 
therefore, the additional 3 per cent associated with the un- 
certainty in the measurement of instrumental resolution, to 
be added in quadrature, does not change the total error sig- 
nificantly. As a check, we repeated the fits with templates 
broadened to the resolution changed by a standard devi- 
ation. For galaxy II we obtained the largest difference (7 
km s _1 ), still negligible with respect to the other sources 
of error. The effect on the uncertainty on the resolution of 
EFOSC2 is smaller and totally negligible. 

The spectra of galaxies Dl, E2, F2, G3, 13, L3, 
M3, R3, S3, U3, D4, H4, and 14 include the lines Ca 
H and K. As noted in T99, even though widely used to per- 
form such measurements (Dressier 1979), these lines have 
been reported to induce a slight overestimate of the velocity 
dispersion (Kormendy 1982; see also Kormendy & Illing- 
worth 1982). Kormendy & Illingworth (1982) suggest that 
the problem may be caused by the intrinsic width of the lines 
or by the steepness of the continuum in that spectral region. 
However, Dressier (1984) suggests that they may be best 
suited for the measurement of large velocity dispersions for 
faint objects, such as Dl. For these reasons, we performed 
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Figure 12. Relative difference in velocity dispersion as measured 
with and without Ca H and K as a function of velocity dispersion 
(upper panel) and average S/N per pixel (lower panel). The error 
bars are the quadratic s um of the systematic and random errors, 
as defined in the Section 4.3. The average difference is 1.4 per cent 
with a scatter of 6.6 per cent. The average difference is negligible 
with respect to the errors (see Table and consistent with the 
scatter found. 



the kinematic fit also by masking the region of Ca H and 
K. In Figure |l^ the variation of a as a function of veloc- 
ity dispersion and S/N is plotted. The velocity dispersion 
obtained with Ca H and K, as listed in Table M, is on aver- 
age 1.4 per cent higher than the one obtained by masking 
this region, with a 6.6 per cent scatter. The effect is within 
the estimated error and the scatter is likely to be induced 
by the loss of information caused by masking a significant 
part of the spectrum. We conclude that with this kind of 
resolution, sampling, signal-to-noise ratio, and treatment of 
the continuum, the presence of Ca H and K does not alter 
significantly the measurement of velocity dispersion. 

Similarly, the spectra of Al, C3, T3, and A4 include 
the region of NaD, which may be affected by interstellar 
absorption even in elliptical galaxies (e.g. Dressier 1984). 
We repeated the analysis by excluding the NaD region (see 
Figure |l3"| ) and found an average velocity dispersion higher 
by 2.3 per cent. As above, we conclude that the presence of 
NaD in the spectrum does not alter significantly the result, 
and used in Table |^ the values found by using the entire 
spectral range. 

The total error on the redshift, taking into account the 
error on wavelength calibration, is less than 0.0005, when the 
GHFF fit was performed and less than 0.002 in the cases 
where the redshift has been measured by identifying the 
main spectral features. 



4.4 Aperture correction 

Spatially resolved kinematics measurements in nearby galax- 
ies (e.g., Capaccioli et al. 1993; Carollo & Danziger 1994a,b; 
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Figure 13. Relative difference in velocity dispersion as measured 
with and without NaD as a function of velocity dispersion (upper 
panel) and average S/N per pixel (lower panel). The error bars 
are the quadratic su m of the systematic and random errors, as 
defined in the Section 4.3. The average difference is 2.3 per cent, 
negligible with respect to the errors (see Table Q). 



Bertin et al. 1994) have shown that the velocity dispersion 
varies (generally it declines) with radius. For this reason, 
central velocity dispersion is generally different from the 
velocity dispersion measured from the spectrum integrated 
over the entire galaxy, as is generally available at interme- 
diate redshift. A correction is thus required. As in T99, we 
model the velocity dispersion profile with a power law, 



(11) 



with —0.1 < d < 0. The desired correction can be computed 
numerically from 



a 2 (A) 



2irrdrcr (r)DV(r) 



(12) 



Here we have assumed that the value obtained by measuring 
a 2 within an aperture A is the average of <r 2 (r) weighted 
by the luminosity density, modelled as an r 1//4 law appro- 
priately normalized and indicated by DV(r). With respect 
to T99, the way we compute the correction has been im- 
proved by taking into account the effect of seeing. The ef- 
fect of seeing is to smear the dependence of the correction on 
re, whenever the two quantities are comparable. We com- 
puted the correcting factor B(d) = cr/fg^ff 1° an equiva- 
lent aperture of radius re/8 (see j0rgensen et al. 1995; T99) 
for a range of values of seeing (0!.'8 — l'.'2), effective radius 
((y/5 — 2"), and number of lines used. The average correct- 
ing factor, computed as B = [B(— 0.1)+jB(0)]/2, ranges from 
1.055 to 1.075. The uncertainty related to the intrinsic va- 
riety of velocity dispersion profiles (modelled by d) is larger 
than this interval. Therefore, we adopt the mean correcting 
factor B = 1.065 ± 0.037, where the error is estimated as 
B = [B(-0.1) - B(0)]/2V3 as in T99. The values obtained 



for the central velocity dispersion are shown in Figures 
and H together with the spectrum of each galaxy. 



4.5 Relative flux calibration 

The instrumental response was measured at the beginning 
and at the end of every night, by observing spectropho- 
tometric standard stars at parallactic angle through a 5" 
slit. The response was stable from night to night within 2-3 
per cent over the entire spectral range when the same star 
was used for calibration. Different stars provided similar re- 
sponses within 5 per cent. We used the average response 
function to correct the spectra. We estimate the relative 
flux calibration uncertainty to be ~ 5 per cent, dominated 
by the systematic differences found when different stars are 
used. 



5 SUMMARY 

In this paper we have presented the surface photometry and 
spectroscopy of a sample of HST selected early-type galax- 
ies. The sample is chosen mainly based on morphology and 
colour. As described in Section H, our sample is a fair sam- 
ple of the early-type galaxies found in the Medium Deep 
Survey, with respect to morphological selection, colour dis- 
tribution, bulge-to-total luminosity distribution, and lumi- 
nosity distribution. As already noticed by other authors (see, 
e.g., Totani & Yoshii 1998), at any given redshift early- type 
galaxies span a significant range of colours, from the red 
envelope of old cluster ellipticals bluewards (Figure ^ . 

In Section ^| we have presented the surface photome- 
try of 35 objects in the redshift range z — 0.11 — 0.65 (see 
Figure ^). The photometric structural parameters (effec- 
tive radius and effective surface brightness) have been de- 
rived with two different techniques: isophote profile fitting 
(r 1,/4 and r 1 / 4 + exponential disc) and 2 dimensional fitting 
(r 1,/4 only). The value of the structural parameters changes 
significantly (r.m.s scatter of 15 per cent in re) with the dif- 
ferent modelling, but the combination of effective radius and 
surface brightness that enters the FP is remarkably stable 
(r.m.s. scatter 0.03). 

In Section |i| we have described spectroscopic measure- 
ments obtained at the ESO-3.6m telescope. Redshifts are 
obtained for all the 35 objects with measured surface pho- 
tometry. Using Montecarlo simulations, we have studied ex- 
tensively the systematic errors in the measurement of a, due 
to instrumental resolution and finite S/N. The systematic 
errors are found to be significant for values of a comparable 
to the instrumental resolution and for small S/N. Based on 
these simulations, we measure velocity dispersion only on 
the 28 spectra with average S/N per pixel larger than 12, 
and we reject all values of a < 125 km s~ 4 . In the end, we 
have obtained robust determination of velocity dispersion 
for 22 galaxies. 

A range of stellar templates is used to recover the inter- 
nal velocity dispersion; the best-fitting template is found to 
span the entire range G4-K0 III. In this way, we are able to 
investigate the systematic errors induced by the mismatch 
between the galactic stellar population and the stellar tem- 
plate used. The average error on internal velocity dispersion 
is found to be 8 per cent (plus 5 per cent systematic). No 
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Figure 14. Distribution of redshift for the sample galaxies 
(empty histogram). The galaxies with measurable velocity disper- 
sion (S/N > 12) are plotted as hatched histogram. The velocity 
dispersions smaller than our resolution li mit were rejected based 
on Montecarlo simulations (see Section The galaxies with 

measured velocity dispersion are plotted as filled histogram. 



significant offset in the average velocity dispersion is found 
when Ca H and K and NaD lines are masked in the kine- 
matic fit. 
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APPENDIX A: 
GALAXIES 



IMAGES OF NOTEWORTHY 



In this appendix we show the images of the two noteworthy 
galaxies E3 and F3. Galaxy E3 shows a clear spiral pattern 
in the residuals. The center of galaxy F3 is obscured by a 
dust lane; the luminosity profile has been measure by cor- 
recting the dust extinction as described in Section H. The 
images of all the gala xies are available at the MPS world 
wide web site at URL http://ar chive. stsci.edu/mds/| . 
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Figure Al. WFPC2 images and 2D fit residuals of secondary target F3. The horizontal bar is 1 arcsec long. 
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Figure A2. WFPC2 images and 2D fit residuals of galaxy F3. The first line shows the F606W images, the image corrected for dust 
extinction (see Section H), and the residuals from the 2D fit. The same images are shown for filter F814W on the second line. The 
horizontal bar is 1 arcsec long. The noise in the central part of the images on the second and third column is higher than the rest, 
because of the additional noise due to the internal extinction correction. 
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